Itinerant and localized magnetism on the triangular lattice: 
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We study the interplay between correlation, itinerant ferromagnetism and local moment formation 
on the electron doped triangular lattice of sodium cobaltates Naa;Co02. We find that strong cor- 
relation renormalizes the Stoner criterion and stabilizes the paramagnetic state for x < — 0.67. 
For X > Xc, ferromagnetic (FM) order emerges. The enhanced Na dopant potential fluctuations 
play a crucial role in the sodium rich phases and lead to an inhomogeneous FM state, exhibiting 
nonmagnetic Co^^ patches, antiferromagnetic (AF) correlated regions, and FM clusters with AF 
domains. Hole doping the band insulator at x=l leads to the formation of local moments near the 
Na vacancies and AF correlated magnetic clusters. We explain recent observations by neutron, fj.SK, 
and NMR experiments on the evolution of the magnetic properties in the sodium rich phases. 

PACS numbers: 74.25.Ha, 71.27.+a, 75.20.Hr 



Sodium cobaltates, Naa:Co02, are electron doped tran- 
sition metal oxides with a layered, hexagonal lattice 
structure. The Co^+ has a 3d^ configuration with 5 d- 
electrons occupying the three lower t2g orbitals. At a Na 
density x, the average Co valence is Co^^~^-*, evolving 
from an open shell Co'^^ with low-spin state S = 1/2 
at a; = to a closed shell Co'^"'" with a low-spin state 
S ^ a,t X — 1. The cobaltates provide an almost ideal 
low-spin lattice fermion system for studying the physics 
of correlation and geometrical frustration in the charge, 
spin, and orbital sectors. 

Since the discovery of unconventional superconductiv- 
ity near x — 0.3 upon water-intercalation in this ma- 
terial pj], a broad spectrum of experiments have been 
performed, yielding a rich and complex phase diagram 
with many unexpected and novel properties. Magnetism 
plays an essential role in the sodium rich region with 
X > 0.5. The metallic transport coexists with a Curie- 
Weiss magnetic susceptibility leading to a novel "Curie- 
Weiss metal" [3]. Neutron scattering and NMR exper- 
iments find strong in-plane ferromagnetic (FM) fluctu- 
ations at high temperatures and magnetically ordered 
states emerge for x > 0.75 [1]. At a; = 0.82, neu- 
tron scattering determined that in-plane FM order and 
inter-layer antiferromagnetic (AF) order develops below 
a 3D Neel temperature T/v — 20K There is increas- 
ing experimental evidence for unexpected strong correla- 
tion effects as the cobaltates approach the band insulting 
limit at X — 1. Recent /iSR experiments discovered the 
emergence of AF correlated magnetic clusters and local- 
ized magnetic moments for x > 0.96 [B]. Interestingly, 
the high thermoelectric power of the sodium cobaltate is 
found in the sodium rich region which is likely to have a 
related magnetic origin [3| . The nature of the magnetism 
in the sodium rich region has been the focus of several 
theoretical work 0,ll,M[i3|- 

In this paper, we study the magnetic properties, both 
itinerant and localized, of strongly correlated electrons on 



the triangular lattice, and present a theoretical descrip- 
tion of the novel magnetism observed by experiments 
in the sodium rich phases. Strong correlation plays an 
important role in the electronic structure of the cobal- 
tates. It has been shown recently [ill [l^ . [isl [3] that 
correlation-induced corrections to the crystal field split- 
ting and bandwidths of the LDA calculations [l5l| can 
produce a single renormalized band of mostly aig char- 
acter at the Fermi level as observed by angle resolved 
photoemission (ARPES) [H, 0, [H, Here we show 
that the description of the magnetism in the cobaltates, 
too, requires a proper account of the large Coulomb re- 
pulsion U at the cobaltate sites. Indeed, in the LSDA-I-U 
theory [23|, which accounts for the correlation effects in 
the Hartree-Fock approach, the LDA paramagnetic (PM) 
state is always unstable and the ground state is a fully 
polarized ferromagnet at all doping for U as small as 
2eV. This disagrees with experiments that the in-plane 
FM order does not emerge until x > 0.75 and that a sin- 
gle spin-degenerate Fermi surface of the Luttinger area 
is observed by ARPES for x < 0.72 We show that 
the weak-coupling Stoner instability is unphysical when 
correlation is strong. Using a strong-coupling Gutzwiller 
projection approach, we show that the Stoner criterion is 
strongly renormalized and the spin dependent self-energy 
scales with the average kinetic energy instead of U. The 
PM phase is in fact stable against itinerant FM below a 
critical electron doping Xc — 0.67 above which in-plane 
FM order emerges. 

In the sodium-rich region, we show that the enhanced 
electrostatic potential fluctuations due to the disordered 
Na dopants lead to the coexistence of localized and itin- 
erant electronic states with inhomogeneous FM order, 
exhibiting nonmagnetic Co^+ patches, AF correlated lo- 
cal regions, and FM clusters with AF domain walls. For 
very high Na doping, the dilute Na vacancies enhance 
the strong correlation effects by increasing the localiza- 
tion tendency of the carriers |21| . We consider the case of 
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a few holes (Na vacancies) doped into the band insulator 
at X = 1, and find that the hidden correlation effects are 
brought out upon the slightest amount of doping. Specif- 
ically, a single hole/Na vacancy induces an S* = 1/2 local 
moment. We address the interactions between the local 
moments and the evolution into magnetic clusters and 
eventually to the macroscopic FM ordered state with de- 
creasing Na doping and provide a theoretical description 
of neutron [l], fj.SR ^ and NMR experiments [22|- 

We start with the Hubbard model for the relevant low- 
energy quasiparticle band of approximate aig character, 
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E 



(1) 



where creates a hole and U (~ 5eV for the cobaltates 
[l^ ) is the on-site Coulomb repulsion. To model the aig 
band, we consider up to 3rd nearest-neighbor hopping 
(ii,t2,i3)=(-0.202, 0.035, 0.029)eV The hole den- 



sity Ui 



\^Xi with Xi the electron doping, is 



fixed by the chemical potential /i. Due to the small direct 
Co-Co overlap and the large U and the 90° 0-Co-O bond 
angle, the AF superexchange J in the cobaltates is small 
[24| , consistent with the value J ^ 5meV determined by 
inelastic neutron scattering We thus focus on the 
in-plane magnetism of the kinetic origin. To study the 
3D magnetic ordering transition at finite temperature, a 
small inter layer exchange coupling is needed. 

In the weak-coupling HF theory, itinerant FM is due 
to the Stoner instability, i.e. the divergence of the uni- 
form susceptibility x = Xo/(l — UNp/2), where xo is 
the free-electron value and Np is the DOS at the Fermi 
level. The large DOS of the cobaltates would lead to FM 
order for all x for U greater than a value less than 2eV 
[20 |. clearly inconsistent with experiments. The failure 
of the HF or LSDA-I-U theory lies in the spin-dependent 
self-energy correction that scales with U . This is unphys- 
ical for U larger than the bandwidth. For example, for 
large ?7, the system can simply avoid paying the energy 
penalty for double occupation in Eq. ([T]) by reverting to 
a fully spin-polarized state that involves only the kinetic 
energy. This is the physics behind the Nagaoka theo- 
rem: on square lattices, the ground state of the infinite-U 
Hubbard model doped with a single hole is a fully po- 
larized FM. At finite hole density, the Nagaoka state is 
lower in energy than the Gutzwiller projected PM state 
at low doping [2a | . Interestingly, on the triangular lattice 
the kinetic energy is frustrated in the sense that hopping 
around an elemental triangle picks up a negative sign, the 
Nag aoka state is not the ground state for a single hole 
27| . To study FM at finite electron doping, we consider 
the large-U limit captured by the projection of double 
occupation. To make analytical progress, we treat the 
latter by Gutzwiller approximation (GA) in the grand 
canonical ensemble where the projected wave function 
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FIG. 1: Renormalized a\g quasiparticle dispersion (a) and 
the DOS (b) in the PM state at x=0.3. (c) The uniform 
susceptibility showing the renormalized Stoner instability at 
Xc — 0.67. (d) The magnetization phase diagram. 



can be written as 



(2) 



Here is an unprojected Slater determinant state and 
Uia is a spin-dependent local fugacity that maintains the 
equilibrium and the local densities upon projection. The 
GA is then equivalent to minimizing the energy of the 
renormalized Hamiltonian, 

Hga = ^ 9tjUjc\„Cja + ^ ei„{c\^Ci„ - n,„), (3) 

i,j,(T 1,(7 

where the Gutzwiller renormalization factor, 



(vI/o|cLc,.|*o) 



(1 - njo-)(l - nj„y 
and is determined by d{HQA) / dni„ — 0, 



(4) 



_ „ E S^ijUj {c\„Cj^ + h.c.) - Em. (5) 



Physically, Si^ is the local kinetic energy per doped spin- 
cr electron measured relative to the average over both 
spins, = - I]a(l ^ nic,)eicy. 

In the uniform phase, riia- = n^, gfj = g"' , Eia = s^, 



1 



(1 - n,) 



^EkafiEka + £a ~ ^J-) - eo, (6) 



where / is the Fermi function and Eka = .9'^[2ti(cos ky + 
2 cos cos ^) + 2t2(cos VSk^ + 2 cos cos^) + 
2^3 (cos COS \/3kx cos ky)] is the renormalized quasi- 

particle dispersion. The quasiparticle band and the DOS 
in the PM phase are shown in Fig. la at j; = 0.3. 
The bandwidth reduction is due to the Gutzwiller fac- 
tor g°' < 1. The peak in the DOS near the band top is a 
property of the cobaltate aig band captured with up to 
third n.n. hopping. 
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FIG. 2: (Color online) The spin (black and white arrows) 
and charge (local doping Xi) distribution in a typical Na(l) 
and Na(2) realization at a; = 0.8 (a) and x — 0.97. (c) The 
one-particle spectral intensity (top) and its spin resolved com- 
ponents at X = 0.8. 

To drive a renormalized Stoner theory for the instabil- 
ity of the PM phase against uniform FM order, we cal- 
culate the uniform magnetic susceptibility x — dM/dh, 
where M = (rif — n|)/iB is the magnetization due to 
an infinitesimal external magnetic field h. Including the 
leading /i-dependence of Eg- in Eq. ([6]) we obtain 



X 



Xo 



1 - KNf ' 



K 



{Ek^ + e + £o) , (7) 



where Np is the renormalized DOS and K measures 
the energy per electron at the Fermi level. The latter 
plays the role of the effective interaction strength, replac- 
ing the Hubbard-U in the Stoner susceptibility. Fig. lb 
shows the doping dependence of x which diverges at 
Xc — 0.67, corresponding to the renormalized Stoner cri- 
terion KNp — 1, beyond which FM order develops at 
T = 0. The self-consistently determined magnetization 
shown in Fig. Ic indicates a sharp transition into the fully 
polarized FM, half-metallic state. The presence of the 
DOS peak near the band top is important for the emer- 
gence of the in-plane FM order. For example, with only 
nearest neighbor hopping, the FM order is completely 
suppressed by strong correlation 0]. We obtained the 
same results from the three-band Hubbard model of the 
t2g complex. The finding of the fully polarized FM state 
at large x is consistent with the large FM moment of 
about 0.13/i_B per Co site a.t x — 0.82 observed by neu- 
tron scattering [4] where 0.18/iB corresponds to a fully 
polarized FM state. 

Now we turn to the spatially unrestricted solutions of 
Eqs. (I3][6]) and study the effects of the Na dopants. It 
is known that Na orders at a; = 0.5 into y/S x 2 su- 
perstructure (2§|, which has been shown [23| to play an 
important role in alleviating geometrical frustration for 
the emergence of the 2x2 AF ordered state observed 



by neutron scattering f25|. For x ^ 0.5, Na ions are dis- 
ordered, and occupy randomly the preferred Na(l) and 
Na(2) sites directly above/below the Co or the middle 
of the Co triangle respectively, with a ratio of 1 : 7 at 
X = 0.8 [ll]. The electrostatic potential is described by 
adding to the Hamiltonian dT]), 



1 



, (8) 



where Vd is the dopant potential strength, dz — a the set- 
back distance of Na to the Co plane, and V the long-range 
Coulomb interaction that must be included to account for 
carrier screening. Figs. 2a and 2b show the typical charge 
and spin distribution of the inhomogeneous FM state at 
X — 0.8 and x = 0.97 on triangular lattices of 32 x 20 
sites, with V — 0.2eV and Vd — 0.6eV. Note that at such 
strength of {V,Vd), the ordered Na at a; = 0.5 only in- 
duces a weak charge modulation ^sf. In contrast, the 
random distribution of Na at x = 0.8 leads to large fluc- 
tuations in the local electrostatic potential, causing the 
localization of the electrons and the formation of the non- 
magnetic Co'^"'" clusters where Xi ~ \. This result is in 
line with NMR: the presence of Co'^+ for x > 0.65 but not 
at X = 0.5 [13, [2^. The one-particle spectral intensity in 
Fig. 2c clearly demonstrates the coexistence of localized 
and itinerant band-like states. Interestingly, AF corre- 
lated regions emerge ai x = 0.8 in locally underdoped re- 
gions in Fig. 2a. In these regions, x is much smaller and 
the kinetic AF correlation imbedded in the Gutzwiller 
factor in Eq. (3) prevails due to the alleviated AF frus- 
tration by charge inhomogeneity, as proposed for the "0.5 
phase" 23|. As a result, the average magnetic moment 



in Fig. 2a is reduced from the fully polarized value to 
about 0.13/iB per Co site, in qualitative agreement with 
the finding of neutron scattering [3]. For stronger po- 
tential fluctuations, FM clusters with AF domain walls 
emerge, as seen more prominently at a; = 0.97 in Fig. 2b. 
These glassy behaviors arise because the localized states 
formed out of the majority spin band are pushed above 
Ep and occupied by holes, as shown in the spin-resolved 
spectral intensity in Figs. 2c. 

To further study the localized magnetism, we consider 
the limit x = I and ask what happens when a few holes 
are doped into the band insulating NaCo02 by the Na 
vacancy. The hidden correlation effects are brought out 
by the slightest amount of doping. Fig. 3a displays the 
case of a single hole added by a Na(2) vacancy. Instead 
of adding the hole into the minority of the spin-polarized 
bands, localized states are created and pinned near the 
Fermi level (Fig. 3a) to accommodate the doped hole, 
leading to the formation of the spin- 1/2 local moment 
distributed near the Na vacancy as shown in Fig. 3b 
for Na(l) and Na(2) respectively. The localized states 
induced by the Na electrostatic potential |21J are spin- 
split by Ai = - e,x = {l/2xi) J^j 9lj{ti]c\^Cj^ + h.c). 
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FIG. 3: (Color online) (a)Spectral intensity of one-hole doped 
by a Na(2) vacancy at Vd = l.OeV. (b) Spin and charge dis- 
tribution of the S=l/2 local moment near a Na(l) and a 
Na(2) vacancy, (c) The local moment "binding energy" Et 
as a function of Vd induced by a Na(l) or Na(2) vacancy at 
Vd = 0.6eV. (d) Spin/charge distribution around two Na(l) 
and two Na(2) vacancies showing AF correlations between the 
local moments. 

The latter has a localized profile whose amplitude Ei, = 
Amax — Amin is uscd as a measure of the binding energy of 
the local moment. In Fig. 3c, E/, is plotted as a function 
of the bare Na potential Vd- The local moment develops 
{Eb > 0) for Yd as small as 0.2eV. Because the Na(l) va- 
cancy is directly above a Co site, Eb is enhanced, making 
it easier for the local moment to form near Na(l) vacan- 
cies. For even smaller Yd, it becomes difficult for our 
finite size numerical calculations to discern the localized 
states and the values of E^, leaving open the possibil- 
ity of self-trapped, spontaneous local moment formation 
without the Na potential. 

When the Na vacancies are isolated, the spin-1/2 local 
moments behave as free moments, contributing to sig- 
nificant spin entropy in the Na-rich part of the phase 
diagram. As the vacancy density increases, i.e. as the 
average x reduces, the local moments begin to overlap 
and their interactions become important. We find that 
nearby local moments have AF correlations. As shown in 
Fig. 3d, since the AF frustration is alleviated by the asso- 
ciated charge inhomogeneity, two Na(2) or Na(l) vacan- 
cies induce two AF correlated local moments with zero 
net magnetization. As the local density of Na vacancy 
increases, FM clusters develop which eventually evolve 
into the macroscopic FM state. Our findings provide 
a complimentary description of the evolution from local 
magnetic clusters to macroscopic FM state observed by 
recent /xSR experiments [1| and interpreted in terms of 
a change in the Co'^"'" smn state that involves the higher 



Co-3d Eg orbitals 

We conclude with a discussion of an outstanding puzzle 
in the ARPES experiments. A large hole- like FS is ob- 
served at a; = 0.8 with a volume much larger than what is 
expected by the Luttinger counting in a PM state ^M^. 
A natural explanation is that this corresponds to the mi- 
nority spin band of the itinerant half-metallic state with 



in-plane FM order. Such an interpretation, although in 
line with that of the neutron scattering experiments 0] , 
implies a filled majority band below the Fermi level which 
has yet to be detected by ARPES. Our study suggests an- 
other possible scenario where the in-plane FM order asso- 
ciated with the 3D A-type AF order does not develop on 
the 2D surface at finite temperatures. The larger than 
expected FS comes instead from the loss of the doped 
electrons due to Na-induced localization and are "taken 
out" of the t2g bands. For example, a substantial frac- 
tion of the doped electrons on the surface can be localized 
to form the nonmagnetic Co'^"'' similar to Fig. 2a. This 
picture is appealing given the recent finding by NMR 
of valence disproportionation associated with significant 
Co'^^ formation for 0.65 < a; < 0.75 [l^]. Interestingly, 
the same picture within the present theory suggests that 
in the extremely sodium rich phases x > 0.96, it is the 
localization of a fraction of holes near the Na vacancies 
that gives rise to the local magnetic clusters observed by 
/iSR experiments [^. Further studies are clearly needed 
to better understand the coexistence of itinerant and lo- 
calized magnetism in the sodium rich cobaltates. 
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